Implications of su(2) symmetry on the dynamics of population difference in the 

two-component atomic vapor 
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We present an exact many body solution for the dynamics of the population difference N2—N1 in- 
duced by an rf-field in the two-component atomic cloud characterized by equal scattering lengths. 
' We show that no intrinsic decoherence occurs for N2 — Ni , provided the exact intrinsic su(2) sym- 

metry holds. Decoherence for N2—N1 arises when deviations from the symmetry exist either in the 
scattering lengths or in the confining potentials. We suggest testing the evolution of N2 — Ni by 
conducting echo-type experiments. 
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Successful achievement of Bose-Einstein condensation in atomic traps M stimulated strong interest in the mech- 
anisms of dissipation and decoherence in the confined atomic clouds (see in 0). Many intriguing questions are 
associated with the finiteness of the number of bosons N forming the Bose-Einstein condensate (BEC) in the trap. 
It has been realized that the global phase may experience collapses and revivals as long as N is finite 

The experimental study of the temporal correlations of the relative phase of the two-component 87 Rb BEC has been 
conducted by the JILA group [g by means of measuring the dynamics of the population difference N 2 — Ni between 
the components. No decay of the global relative phase has been detected on the time scale of the experiment < 100ms. 
The question was then posed in |8|] as to why the phase correlations are so robust despite an apparent fast relaxation 
q . of other degrees of freedom such as, e.g., the relative motion of the two condensates. 

In the present work, we address the question |8| of the relative phase correlations. We show that if the two- 
component system obeys the intrinsic su(2) symmetry, as is almost the case for 87 Rb ||, no intrinsic decoherence 
■ exists for the operators which are the generators of this symmetry. In particular, the population difference falls into 
this category. This result is general and does not depend on a particular form of the Hamiltonian. On the contrary, 
QO ' the- nature of the decoherence arising due to deviations from the symmetry is not universal, and it depends on details 
I/"") \ of the Hamiltonian. We discuss how the nature of the decoherence of N 2 — N\ can be tested in the echo-type 
experiment (see, e.g., in JlfH). Specifically, we suggest modifying the JILA experiment j^] by introducing an extra 
pulse between the two n/2 pulses employed in ||. The echo should then be seen by the readout pulse (the last 
one) at the time t e = 2r where r is the time delay between the first two pulses. We note that the atomic echo in 
optical lattices, predicted by one of the present authors has been observed very recently Q under, however, the 
conditions when the inter-particle interaction is not important. 

We consider the two-component bosonic system with no exchange interaction between the components (no intrinsic 
inter-conversion between the species). This system can be described by the second quantized Hamiltonian H [ p"3[ . 
We note that this Hamiltonian has been employed for studying the dynamics of N 2 — Ni in 87 Rb in the mean 
field approximation (MF) pjjl5]|. In what follows we will obtain the class of exact solutions of the full many-body 
quantum mechanical problem in the case of the su(2) symmetry. We represent H [l5| ] as 

H = H s + H A + e J z +Hp, (1) 
H s = y"dx{*Ii/ ls * 1 +^ i/l ^ 2 + |(^t^ i + ^t^ 2)2 + | ( ^t^ i _^t^ 2) 2| j (2) 

Ha = J d^{^\H la ^ 1 -^lH la ^ 2 + 9 -[{^ 1 ) 2 ~{^ 2 ) 2 ]}, (3) 

Here $ \ is the second quantized Bose field of the z-the component (i = 1,2); H$ is symmetric and Ha is anti- 
symmetric with respect to exchanging the components $ li2 — > ^2,1; Hp describes the effect of the rf-field ]|,|ld],|l5]] ; 
f-o = const is the detuning of the external field hCl(t) (taken in the rotating wave approximation ) from resonance 
between the components; fl(t) stands for the corresponding Rabi frequency treated as an envelope of the rf-pulses 
§||; H u = (Hi + H 2 )/2, H la = (Hi - H 2 )/2 ; Hj = -h 2 S7 2 /2m + U^x) (j = 1,2), stands for the one-particle 
part which includes the kinetic energy and the trapping potential; the binary collision terms are taken in the contact 
form, where 
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7T?i 2 , 7T?l 2 . 2'KTl' 

.g = (ai + a2 + 2ai 2 ), g = (oi + a 2 - 2ai 2 ), .g = (01-02) (5) 

mm m 

are related to the corresponding scattering lengths a\, a 2 , eti 2 . We have also introduced the following operators 

J z = i /" dx(¥j* 2 - *l*i); J+ = J dx$ys x \ J- = [dx&l*2. (6) 

It is important that the operators (^|) represent the su(2) algebra of the angular momentum operators. Indeed, it is 
easy to see that these operators obey the standard su(2) commutation relations 

[J Z ,J+] = J+, [J Z ,J_]=-J_, [J+, J_] =2J Z , (7) 

provided the field operators obey the Bose commutation rule [^(x), ^j(x')] = SijS(x — x'). The operator J z = 
(N 2 — Ni)/2, with N12 standing for the total number of bosons in the first and the second components, respectively, 
can be viewed as the z-component of the angular momentum operator; J x = (J+ + J_)/2, J y — (J+ — J_)/2i are 
the x,y-components, respectively. We note that the formalism of the effective angular momentum operators has been 
successfully applied in many areas including optics (see, e.g., in Refs. J|jl^]). 

In the work |J it has been emphasized that the two component 8 Tlb realizes the su(2) intrinsic symmetry to a 
very good approximation due to almost perfectly satisfying the relations a\ = a 2 = ai 2 as well as U\ — C/ 2 . Thus, 
g = g' = and H± a = in Eqs.(0), (g), (||). Then, the part Ha @ is zero, and in the case e = 0, Hp = , 
the total Hamiltonian is just H$ (0) which commutes with all the su(2) operators (|^) even though g ^ 0. This 
situation corresponds to an exact su(2) symmetry Q in the trap |L7| . 

Let us, first, consider the effect of violation of the su(2) symmetry by the external rf-field (no intrinsic violation). 
So we choose g' = 0, 5 = 0, Hi a = in the Hamiltonian (0)-(Q). It is important to emphasize that in this case, 
the full Hamiltonian (]l|) forms a closed algebra with the operators (||). Accordingly, the Heisenberg equations for 
J+,-, z are 

ih j z = m W j + _ j_ t ihj + = -e J++htt(t)J z , ihj- = e J--hSl*(t)J z , (8) 

where the su(2) commutation relations (Q) have been employed. We note that these equations are exact and linear 
in ■ Accordingly, there is no dissipation for these quantities ( unless the parameters 60 an d £l(t) fluctuate 

due to external noise). 

The above conclusion of the absence of the intrinsic decoherence for the population difference in the case of intrinsic 
su(2) symmetry is the main result of this work. It is important that it is not limited by the specific form of the 
Hamiltonian In fact, any form of H where U\ — C/ 2 and the interaction part can be expressed solely in 

terms of the total density operator p = W-^i + V I / 2 V I' 2 will lead to the same conclusion. This is so because H and 
J z t ± form a closed algebra in this general case as well. 

Now let us consider particular cases. First, we study the effect of the two short pulses investigated experimentally 
in Ref. In the initial state \t — 0), all N = Ni + A 2 bosons occupy only one component (e.g.,i = 2), so 
that (t = 0|J z |f = 0) = N/2 ||. The effect of a short pulse can be described by the unitary transformation 
U = cxp(— i J dtHp) in the sudden approximation (from now on we choose units h = 1). Taking into account Eq.(|]), 
we find 

U = e -i(v'J„+v"Jy) > (9) 

where v' — jdtfl'(t), v" = JdtQ"(t), with Q'(t) and Q"(t) being the real and imaginary parts of fi(i) in Eq. ([!]). 
As a result of the pulse taking place at some time moment t, any operator 0(t) before the pulse transforms into 
U'0(i)U after the pulse. Specifically, the operator J z changes as 

J z -> U f J z U = cosvJ z + — (v'J v - v"J x ) , (10) 

v 

where v = V v' 2 + v" 2 . This relation is a consequence of @j. Let us choose v = \J v' 2 + v" 2 = tt/2, which corresponds 
to the 7r/2-pulse If the first 7r/2-pulse was imposed at time t — 0, and the — 7r/2-pulse was produced at t > 0, 
the state just after the second pulse can be represented as \t) = W exp(—iHt)U\t = 0) = exp(— iWHUt)\t = 0). 
Then, the population difference N2 — N\ = 2J Z after the second pulse is 

2(t\J z \t) = 2(t = 0\e ieotuU * u J z e- ieotu1j * u \t = 0) = N cos(e t) . (11) 
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where the explicit form (0)-(|4|) (for H with g' = 0, g = 0, H\ a = 0) as well as the conditions (t — 0| J x ,y\t = 0) = 0, 
and E g. (]10| ) have been employed for v = ir/2. We also took into account that [Jj, Hg] = 0. 

Eq. represents an exact many body solution in the presence of the binary collisions for the case of exact intrinsic 
su(2) symmetry. The collisions, causing dissipation of the normal modes, do not affect the dynamics of the su(2) 
generators Jj and their observables, as long as the intrinsic su(2) symmetry holds M|. As a matter of fact, in the 
above analysis we never assumed the presence of the Bose-Einstein condensate. Thus, the solution ( |ll| ) is formally 
true in the normal phase as well. In some sense, the intrinsic su(2) symmetry protects the population difference 
from decoherence. This property can be employed in quantum computations ]19|] , where decoherence is usually the 
main fundamental obstacle. 

The goal of the present paper is elucidating the role of the intrinsic su{2) symmetry rather than discussing all 
possible mechanisms of the decoherence, which arise due to violation of the symmetry. This will be done in a separate 
publication. Nevertheless, let us discuss some most obvious sources of the decoherence. First, we comment on the role 
of instrumental noise which is always present. Referring to the solution (|TT|) , the noise may produce fluctuations of eo . 
In general, these can be of two types: i) temporal €i r (t) with some short correlation time ; ii) ensemble fluctuations 
e r from realization to realization (in the destructive measurements fslJlif ) , so that e r = const during each realization. 
Thus, one can represent eo = eo + &ir(£) + e r , where eo = (eo) stands for the mean over the fluctuations (which 
are assumed to be gaussian). The type i) produces irreversible decoherence characterized by the exponential decay 
(t\J z \t) ~ exp(— Tt) , with r = J Q dt(ei r (t)ei r {0)) « Ti r (ei r (0)ei r (0)} at t > Ti r . The type ii) induces reversible 
decoherence (dephasing) characterized by the gaussian time dependence (t\J z \t) ~ exp(— i 2 /rj) with the dephasing 
rate t^ 1 = (e 2 .) /2 . The reversible dephasing will dominate the evolution if the condition t7 1 <C is fulfilled, 
provided t,} ^> r ir . In other words, the irreversible temporal noise (type i)) is irrelevant if \/ jej) ^S> Ti r (ej r (0)) . Thus, 



the decay is reversible at times < t < r _1 . This can be tested by the echo method [HO 1 2 1 . In order to observe 
the echo, three pulses are required, so that the first 7r/2-pulse takes place at t = 0, then some extra (time-reversal) 
pulse is imposed at t — r > , and, finally, the readout 7r/2-pulse takes place at some time t > t. 

Below we present an exact echo solution in the case of the exact intrinsic su(2) symmetry for the case r c i < t < T^ 1 . 
The time-reversal pulse is taken as the — 7r-pulse. Then, the state just after the last pulse is \t) — U exp(—iH(t — 
t))WW exp(-iHr)U\t = 0) = exp(-iUHU*(t - r)) exp(-iW HUr)\t = 0). We take into account (0), and obtain 
the population difference N 2 — Ni — 2(t\ J z \t) (for g = 0, g' = 0, H\ a = ) after the third pulse as 

N 2 -Ni =JVcos(e (t-2r)). (12) 

The solution ([jj]) gives a complete revival of the initial population difference at t — 2t. Averaging over the type ii) 
noise yields 

N 2 - Ni = Ncos(e (t - 2r))e- (t " 2r)2/r ' . (13) 

Should the temporal noise (the type i)) be present, the echo (at t = 2r ) acquires the suppression factor exp(— 2rV) . 
Thus, the echo can be observed if 2iT < 1. 

An intrinsic decoherence of (t\J z \t) occurs when the intrinsic su{2) symmetry is broken. This may be due to 
deviations of g, g' and H\ a from zero fj"7f . Accordingly, the terms describing coupling of the operators (|6|) to the 



normal modes will enter Eqs.(|8|). We note that the effect of H\ a ^ has been studied recently in Ref. |15|. 

Here, we consider the effect of the shot noise produced by uncertainty in the initial depostion of the total N in 
the trap. In this analysis we take into account that g = 0, g' ^ in Eqs.(Q),(||) for the case ||. We also assume that 
the shot noise is large, so that thermal effects can be neglected, and one can set T = 0. Furthermore, we consider no 
asymmetry in the trapping potential ( Hi a = 0). At T — the main contribution to (t\ J z \t) = (N 2 — Ni)/2 comes 
from such many-body eigenstates where almost all N\ :2 belong to the corresponding condensate single particle states 
(/?oi(x), <po2(x). In this case, the operators 'I'i^x) can be represented as *i = ai</?oi(x), and ^2 = a2<^02(x), 
where a± l2 destroy one boson from the condensate states </?oi(x), <^02(x), respectively. Such an ansatz resembles the 
two- mode approximation widely used previously (see, e.g., in Ref. g). Employing it in Eqs.(|)-(|) and (§) yields the 
Hamiltonian (|l|) as 

H = C + (e + hN)J z + b 2 Jl + + W J_, (14) 

where C commutes with (||), and thereby can be omitted as long as the dynamics of J z .± is concerned; the coefficients 
61, b 2 can be expressed in terms of the original parameters entering Eqs. ([!])- (||). In fact, the form (|l4| ) is general as 
long as only the binary collisions are taken into account, and no exchange interaction occurs. It is important that 
61, b 2 vanish if the intrinsic su(2) symmetry holds. The effect of the shot noise is described by the term ~ b^. 
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In this paper we consider the simplest case of the confinement in a box of the volume V, and neglect the phase 
separation effect ||. Then, provided H la — 0, g — 0, we find h = -Stf/V, b 2 — @ in Eq.(0). Now it can be 
seen that the role of £q is played by £q ~ 2g'N/V in the solutions (Jll|) , (12), provided N — const during each 



measurement M. We assume that the shot fluctuations are characterized by some known variance AN. This yields 
the dephasing rate 



_ x _ V2\g'\AN 



V 



(15) 



in Eq. ( jl3| ) . It is important to emphasize that the shot noise does not disrupt the phase coherence in each realization. 
Therefore, the echo effect described by Eq.((l^) will be observed in this case. It is worth also noting that the trap 
losses may introduce some effective temporal noise into Eq. (|l4"|) . It will yield the exponential factor ~ exp(— t/r^) in 
Eq. (|l3|) , where tl is the confinement time. Thus, the trap losses do not disrupt the coherence because in any realistic 
experiment its duration t should be much shorter than tl- 

Let us estimate the rate of the dephasing induced by the shot noise. In 87 Rb, \g'\/g = 0.03; the quantity 
a pa gN/V, which is a typical value of the chemical potential, can be taken as u ps 10 2 — 10 3 s _1 ( % = 1) 0. Thus, 
t^ 1 pa OMfiAN/N pa (4 - 40)AiV/7V s - 1 . For AN pa y/N, we estimate r d pa (0.3 - 3)s from Eq.(|f), where we 
take TV = 10 4 . For N — 10 6 , this estimate increases by the factor of 10. Such long decoherence times should be 
contrasted with the longest relaxation times 10 — 50ms for the normal modes We note, however, that the effect 
of the shot noise may be deliberately enhanced by making, e.g., AN ~ N. In this case the echo effect can be made 
much more pronounced. For example, if AiV = 0.1N for typical number of atoms N pa 10 4 and u sa 10 3 s _1 , we 
find Td pa 30ms. 

The case of T ^ when the intrinsic su(2) symmetry is slightly broken requires special consideration. This 
will be presented in a separate work. Here we only give a qualitative reason as to why the reversible dephasing 
should dominate the evolution of N 2 — N x between the rf-pulses as long as a significant portion of atoms forms the 
two-component condensate. Indeed, if no external bath is present and Ni^ are conserved separately, the dominant 
evolution of the condensate operators ^(x, t) ~ aj(i)y>oi( x ) m & n y given exact many-body eigenstate is simply given 
by the exponents exp(— flit) , where /ij is the chemical potential of the ith component in the given eigenstate. A 
detailed analysis of this has been given in the case of a single component condensate in Ref. pl| . Thus, temperature is 
simply producing the thermal averaging of the exponent exp(i([ii — (^2)t) over the exact eigenenergies. This produces 
the gaussian decay determined by the thermal fluctuations of the relative chemical potential /i =_Ui — M2- I n other 
words, the normal modes of the condensate do not produce irreversible decay of the global phase |2^] . This is in close 
analogy with the case of the single condensate pl| . Thus, the echo effect should persist at T ^ as well, if BEC is 
present pi| . 

In summary, we have pointed out that the dynamics of the population difference in the two-component atomic 
vapor is controlled by the closeness to the su(2) internal symmetry. In the case of the exact symmetry, no intrinsic 
decoherence occurs for the population difference. This result answers the question posed by the JILA experiment 
P|. In the presence of the condensate, slight deviations from the intrinsic symmetry result in reversible dephasing of 
the population difference which can be revealed in the echo experiment. We suggest considering the global relative 
phase of the two-component condensate, which obeys the intrinsic su(2) symmetry and thereby is protected from 
the intrinsic decoherence, as a viable candidate to be employed for quantum computations. 
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